Quantum Accelerator Modes near Higher-Order Resonances. 
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Quantum Accelerator Modes have been experimentally observed, and theoretically explained, in 
the dynamics of kicked cold atoms in the presence of gravity, when the kicking period is close to 
a half- integer multiple of the Talbot time. We generalize the theory to the case when the kicking 
period is sufficiently close to any rational multiple of the Talbot time, and thus predict new rich 
families of experimentally observable Quantum Accelerator Modes. 
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Present-day experimental techniques afford almost per- 
fect control of the state and time evolution of quantum 
systems, and thus allow observation of phenomena, that 
are rooted in subtle aspects of the quantum-classical cor- 
respondence. In particular, effects of mode-locking and 
nonlinear resonance, that are ubiquitous in classical non- 
linear dynamics, could be observed on the quantum level, 
in the form of unexpected quantum stabilization phenom- 
ena; for instance, in nondispersive wave-packet dynamics 
[H, and in the kicked dynamics of cold and ultra-cold 
atoms. In the latter case, techniques originally intro- 
duced by M. Raizen and coworkers have been successfully 
used to produce atom-optical realizations of the Kicked 
Rotor (KR) model 0, which is a famous paradigmatic 
model of Quantum Chaos. A variant of the KR, which 
was realized in Oxford, had the (Cesium) atoms freely 
falling under the effect of gravity between kicks. Discov- 
ery of a new effect followed, which was named Quantum 
Accelerator Modes (QAM) A natural internal time 
scale for the system is set by the so-called Talbot time, 
and whenever the kicking period is close to a half-integer 
multiple of that time, small groups of atoms are observed 
to steadily accelerate away from the bulk of the atomic 
cloud, at a rate and in a direction (upwards, or down- 
wards) which depend on parameter values. A theory for 
this phenomenon introduces a dimensionless parame- 
ter e, which measures the detuning from exact resonance, 
and shows that the nearly resonant quantum dynamics 
may be obtained from quantization of a certain classical 
dynamical system [15|, using e a.s the Planck's constant. 
This dynamical system was termed the e-classical limit 
of the quantum dynamics, and is quite different from the 
system, which is obtained in the classical limit proper 
fi, — > 0. QAM are absent in the latter limit, and arc ac- 
counted for by e-classical phase space structures. Thus, 
they are at once a purely quantal phenomenon, and a 
manifestation of classical nonlinear resonance; indeed, 
their theory is a repertory of classic items of nonlinear 
dynamics, occurring in a purely quantum context. For 



instance, they are associated with Arnol'd tongues in the 
space of parameters, and are hierarchically organized ac- 
cording to number-theoretical rules 0]. Finally, on the 
quantum level, a deep relation to the famous problem of 
Bloch oscillations and Wannier-Stark resonances 0] has 
been exposed 

Existence of QAM somehow related to other rational 
multiples of the Talbot time ("higher order resonances"), 
than just the half-integer ones, is a long-standing ques- 
tion, that lies beyond the reach of the existing theory. 
Some indications in this sense are given by numerical 
simulations, and also by generalizations of heuristic ar- 
guments 0, which were formerly devised Q in order to 
explain the first experimental observations of QAM. 
In this paper we show that QAM indeed exist near reso- 
nances of arbitrary order. This noticeable re-assessment 
of the QAM phenomenon requires a nontrivial reformu- 
lation of the small-e approximation, in order to circum- 
vent the basic difficulty, that no e-classical limit exists in 
the case of higher resonances. We show that, in spite of 
that, families of rays (in the sense of geometrical optics) 
nonetheless exist, that give rise to QAM in the vicinity 
of a KR resonance. Such "accelerator rays" are not tra- 
jectories of a single formally classical system, but rather 
come in families, generated by different classical systems, 
which provide but local (in phase space) approximations 
to the quantum dynamics. This is remindful of the small- 
h asymptotics for the dynamics of particles, in the pres- 
ence of spin-orbit interactions [lo| . This similarity is by 
no means accidental, because the KR dynamics at higher- 
order resonance may be described in terms of spinors 
0, Ell; thus, the present problem naturally fits into a 
more general theoretical framework, and our formal ap- 
proach may find application in the broader context of 
quantum kicked dynamics, in the presence of spin. 
The dynamics of kicked atoms moving in the vertical di- 
rection under the effect of gravity is modeled by the fol- 
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FIG. 1; (Color online) Momentum distributions, in the time 
dependent gauge, after n = 100 kicks, for different values of 
the kicking period near the resonance r = tt ((a) and (b)), 
and near the resonance t = 147r/13 (c). Red color corre- 
sponds to highest probability. The initial state is a mixture 
of 100 plane waves sampled from a gaussian distribution of 
momenta. Vertical dashed lines correspond to the mentioned 
resonant values of t/2ii. Black full lines show the theoreti- 
cal curves (HIJ, with: (a) T = 2,p = 3,j = 1, A2 = , (b) 
r = l,p = 5,j = 1, Ai = and (c) T = l,p = 5,j = 1, Ai = 
207r/13. Parameter values are: k = O.Svr and J? = 0.126r. 



lowing time-dependent Hamiltonian: 
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Units are chosen so that the atomic mass is 1, Planck's 
constant is 1, and the spatial period of the kicks is 2tt. 
The dimcnsionlcss parameters k,T,ri are expressed in 
terms of the physical parameters as follows: k = n/h, 
T = HTG'^/M, 7j = MgT/{hG), where M,T,k are the 
atomic mass, the kicking period, the kick strength, and 
2tt/G is the spatial period of the kicks. X is the posi- 
tion operator (along the vertical direction) and the kick- 
ing potential V{X) = cos(X) in experiments. Hamilto- 
nian (p]) is written in a special, time dependent gauge 
in which the canonical momentum operator is given 
by P + rjt/T. This choice of a gauge makes (Q) in- 
variant under spatial translations by 27r, so, by Bloch 
theory, the quasi-momentum f3 is conserved. With the 
present units, /? is the fractional part of P. The dy- 
namics at fixed (3 are formally those of a rotor with an- 
gular coordinate 9 = Xmod(27r). Let \tpn) denote the 
state of the rotor immediately after the n-th kick; then 
\ipn+i) = Un\ipn), whcrc the unitary operators C/„ are 
given, in the ^-representation, by: 



(2) 



the propagator of the generalized Kicked Rotor. Multi- 
plication of wavefunctions by cxp(im6), (m G Z) 
generates the discrete unitary group of (angular) mo- 
mentum translations. For special values of r and /3 a 
nontrivial subgroup of such translations commutes with 
the KR propagator. This leads to a special dynamical 
behaviour, called KR-resonance @. We define the order 
of a KR resonance as the minimum index of a commut- 
ing subgroup; or, the least positive integer £ such that 
(|3| commutes with multiplication by exp(z^0). KR res- 
onances occur if, and only if [l^ . r is commensurate to 
27r, and the quasi-momentum /? is rational. Indeed, mo- 
mentum translations by multiples of an integer £ leave 
(Q invariant if, and only if, (i) r = 27rp/q with p, q co- 
prime integers, (ii) £ ~ mq for some integer m, and (iii) 
f3 = v /mp+mq/2 mod(l), with v an arbitrary integer. In 
the following we restrict to "primary" resonances, which 
have m = 1 and £ ^ q, and generically denote f3r the res- 
onant values of quasi-momentum. The KR propagator at 
exact resonance is obtained on substituting r = 2'Kp/q, 
(3 = (3r, and 77 = in ([2]). Using Poisson's summation 
formula, it may be written in the form: 
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2^s/q) , (3) 



where 
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Gs=Gs{p.q,f5r) = -E^"" 



-7rip(i+/3,.)^/g l-Kisl j q 



, (4) 
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so that \Gs\ = g"^/^ Now let r = 2pTT/q+e, (3 = Pr + SP; 
and denote 0„ = (5/3 + 77/2 + ijn. We may write 

Un = U,^,cxp{-i^{~id0 + (3r)'^) cxp{-T<f)ndg) . (5) 

Here, and in the following, phase factors only dependent 
on /3 and n are disregarded. Thanks to eqn.([]T]), eqn.([n]) 
may be rewritten in the following form : 



Gsi:{e~2ns/q-TC^,,) , (6) 



s=0 



For r/ — 0, C/„ does not depend on n, and coincides with 



where: 

^{0) = e-'4(-*^''+^'-)V(f) ■ (7) 

If e is granted the formal role of Planck's constant, then 
operator (Q) has the form of a unitary propagator for 
a generalized free rotor [l3|, so quasi-classical methods 
may be used to investigate the small-e regime. We define 
the e-classical momentum operator / = —ied/dO [3]. 
Denoting k = ke, the e-quasiclassical asymptotic regime 
is defined by e — > at constant fc, /. Using the explicit 
form of the integral kernel for ^ [l3|, the transition 
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FIG. 2: (Color online) Phase portraits of maps JTq^' on the 
2-torus, for (a) T = 2,5s = (-l)"+^7r, = -0.395, tt) = 
1.122, (p,j) = (3,1) (r/27r = 0.475, e = -0.157); (b) T = 
1,5^ = 0,fc = 0.032, rr? = 1.253, (p,j) = (5,1) {r/2n = 
0.502, e = 0.013) and (c) T = 1, 5^ = 207r/13, A; = 0.040, r?? = 
1.455, (p, j) = (1, 1) (r/27r = 0.541, e = 0.016). 



amplitude from 6 ~ 9q at time to 9 — 9n after ?7 kicks 
is given by: 

{0n\Un-i---Uo\eo)^{27rtey''/^ ^ G,„...G3„_,x 

(m,s)eO„ 

Ja Jo 

where m and s are vectors with (n — 1) integer compo- 
nents, fin = X {1, ... , and 

71 

SmASo, ■ ■ ■ ,0n) ^J2^-~kV{9t)+ 
t=l 

+ li9t - Ot-1 - 2TTSt/q - 2mtTT - T(btf} . (9) 

Replacing O in ([8]), and using the stationary phase ap- 
proximation in individual terms in the sum on the rhs of 
(I>|, we find that, at small |e|, (P) propagates along rays, 
which satisfy the equations : 

6*4+1 = 9t + It + T(f)t + 27rst/(7 mod 27r 

It+i ^ It - ~kV'{9t+i) , (10) 

or, defining Jt = It + rcpt + l-KStjq, and (5t = 27r(sf+i - 
st)/g, 

Jt+i = Jt +(5t+T?7-fcy'(6lt+i) , 

9t^i = 9t ^Jt mod(27r) . (11) 

For each value of t, (|~T|) defines a map Tt on the cylin- 
der; however, since the choice of the integers si, S2, . . . 
is totally arbitrary whenever g > 1, such maps do not, 
in general, uniquely define a classical dynamical system. 



The St may be removed by changing variables to = q9, 
but this calls into play the function V'i^djq), which is 
not a single- valued function in i? G [0, 27r], except in 
the case when V(9) is a 27r/(7-periodic function; then 
eqs. (fTT|) reduce to a single map, and the theory pro- 
ceeds essentially identical as in the case q = \. In all 
other cases, exponentially many different maps enter the 
game upon iterating eqs. (| |) . and so no e-classical limit 
proper exists. In spite of that, we shall presently show 
how a stability requirement singles out special families 
of rays, which give distinguished contributions in the 
dynamics, ultimately resulting in QAM. In stationary 
phase approximation, each ray (| \ contributes a term 
g-"/2|det(9Jl)|-i/2exp(iS',,™/e-hz$s) in (0), where S^,^ 
is the action Q) computed along the given ray, col- 
lects phases from the Gsj and from Maslov indices, and 
9Jl is the matrix of 2nd derivatives of ([3 with respect 
to the angles 0i, . . . , 9n-\. Stability of a ray is related 
to the behavior of the prefactor |det(97l)|~^/^ as a func- 
tion of "time" n. 2Jl is a tridiagonal Jacobi matrix, with 
off-diagonal elements equal to —1, and diagonal elements 
given by ~kV"{9t) + 2, where 6t are the angles along 
the ray. For a large number n of kicks, most choices of 
sg{1,.. .,(/}" are essentially random. The same may be 
assumed to be true of the diagonal elements of 9Jt, and 
so 971 has a positive Lyapunov exponent, due to Ander- 
son localization. It follows that |det(9Jl)| exponentially 
increases with n (as may be seen, e.g., from the Herbert- 
Jones- Thouless formula [13] )■ Therefore, such rays carry 
exponentially small contributions, and their global effect 
is determined by interference of exponentially many such 
contributions. In contrast, distinguished contributions 
are given by those rays, whose matrices 9Jl have extended 
states, thanks to absence of diagonal disorder. The sim- 
plest such case occurs when the diagonal elements of 971 
are a periodic sequence. This in particular happens when 
6t is a periodic sequence, and rays are in such cases 
related to stable periodic orbits of certain classical dy- 
namical systems, which are constructed as follows. Let 
St+T = St for some T and all t. Then map J^t (| P pe- 
riodically depends on "time" t; so, for each choice of t' 
with < t' < T — 1 , one may introduce a " map over one 
period" ^jJ^ = J-t'+T-i o ■ ■ ■ o Tt', whose iteration de- 
termines rays (fTT|) at every T-th kick after the t'-ih one. 
As this map is 27r-periodic in J, 0, it defines a dynami- 
cal system on the 2-torus. Systems that way constructed 
with different t' are obviously conjugate to each other, so 
the periodic orbits of any of them one-to-one correspond 
to the periodic orbits that are obtained for t' = 0. As a 
result, to each periodic orbit of (on the 2-torus) a 
ray (I ! ■ is associated, which is periodic in position space; 
therefore, its matrix 9Jt has periodic diagonal elements. 
If the orbit has period p, then the corresponding ray (j , i p 
satisfies J^ij^p^T = Jit + 2it) for all integer I, where j is the 
"jumping index" of the periodic orbit. This is equivalent 
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27r(s(;+p)T - siT)/q - T-qTp + 2tt) and 
so, along such a ray, the physical momentum J/e linearly 
increases (or decreases) with average acceleration 

a = e-i{2^j(pT)-i-AT-r77} , (12) 

where At = T^^ X^J^o^ ^s- Finally, stability of such rays, 
as determined by the behavior of dct(9Jt) as a function of 
"time" n, is controlled by the Lyapunov exponent, and so 
is equivalent to dynamical stability of the corresponding 
periodic orbits |17|. 

In summary: whenever V{9) is not 27r/(3'-periodic, no e- 
classical limit exists for the dynamics Q ; QAM may nev- 
ertheless exist, associated with stable "accelerator rays", 
that are associated with the stable periodic orbits of a 
family of maps of the 2-torus. There is one such map for 
each choice of a periodic sequence in {1, . . . , g}^. The 
simplest choice is J* = ; the relevant map (pT|) . and the 
acceleration formula (fT2|) . are then the same as in the 
case 5 = 1 01 ■ In Fig. Q] wc show numerical evidence 
for QAM associated with the resonances at g = 2, p = 1 
((a) and (b)) and at g = 13, p = 7 (c). Here, the kick- 
ing potential is V{0) = cos(0). For the given parameter 
values, three QAM are clearly detected: two around the 
<7 = 2, p = 1 resonance, one near the q = 13, p ~ 7 one. 
They correspond, via eq. (fT2|) . to stable periodic orbits 
of maps T^^^ with T = 1 and with T = 2. The stable 
islands of these orbits arc shown in Fig. Q 
The present theory suggests an unsuspected richness of 
QAM, associated with the dense set of higher-order res- 
onances. If produced with ideal, infinite resolution, fig- 
ures in the style of Fig. P might reveal that QAM are 
essentially ubiquitous; however, some QAM associated 
with resonances of low order q > I should be observable 
already on the present level of experimental resolution. 
Our numerical simulations have exposed a fine texture 
of seemingly QAM-like structures; on the available level 
of precision, however, most of them are so vague, that it 
is impossible to decide to which resonance they belong. 
Those for which this question could be answered were in 
all cases found to correspond to some stable orbits, in 
agreement with the above theory. On the other hand, 
for a few of the periodic orbits we have computed, no 
partner QAM could be detected. This may be due to 
the fact that, at given parameter values, many different 
orbits coexist, which are related to different resonances, 
hence to different values of the pseudo-Planck constant 
e. The hierarchical rules that determine their relative 
"visibility" are not known at this stage. In general, one 
may expect stronger QAM near lower order resonances, 
yet exceptions are not rare, see Fig.l (c). 
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